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Abstract 

This contribution proposes a new family of methods to 
reduce the Peak-to-Average power Ratio (PAR) in 
Discrete Multi-Tone (DMT) and Orthogonal Frequency 
Division Multiplexing (OFDM) systems. This method 
avoids an use of extra Inverse Fast Fourier 
Transformations (IFFTs) as was done in some previously 
published techniques, but instead is based on 
transformations in the time domain of the original 
sequence available at the output of the IFFT processor. 
The improved statistics of peak power in the optimized 
transmit signal are demonstrated by simulation results. 



1. Introduction 

Multi-carrier modulation has recently being used for 
many applications such as high-speed voice band modems 
and digital subscriber line systems (xDSL). These systems 
use DMT over wire media and OFDM for wireless 
communication. 

The amplitude distribution of a DMT signal with 
random input data is approximately Gaussian for large 
number of carriers. Therefore the DMT signal will 
occasionally present very high peaks. These peaks require 
a high dynamic range of the ADC (analog to digital 
converter) and analog front end in absence of any clipping 
or peak reduction technique. This would result in 
inefficient amplifiers (with excessive power dissipation) 
and expensive transceivers. 

Recently, there has been a variety of proposed methods 
to reduce PAR [2]-[6], but none of these methods is able 
to achieve simultaneously a large reduction in PAR with 
low complexity, without performing an extra IFFT 
calculations, without loss in data rate, without increase in 
the transmitter average power and without decrease in the 
minimum distance or equivalently without loss in the 
noise margin 

This paper proposes a new technique based on [1], 
which can achieve all these goals. 



2. The Proposed Method to reduce PAR 

The basic characteristic of DMT is that the data are 
modulated on N QAM subcarriers X k (k = 0, 1, . . . , N - 1 ) 
each of equal bandwidth. This frequency multiplexing can 
^easily be implemented by using an Inverse Fast Fourier 
Transform (IFFT) in the modulator that transforms the 

subcarrier vector X = [x Q , X x , . . . , X N _ x ]' , with 

Xi=X* N _ M for i = 2,...,N/2 and t denotes the 
transpose and * the complex conjugate, into the time 
domain. This results in the discrete time representation 

x = [*o>*i with 

1 &** 

X k e " * = <U...,tf-l (1) 

The sequence xis real and is often referred to as a 
DMT symbol. It's well known that for large N, the output 
time vector distribution can be well approximated by a 
Gaussian distribution and has a relatively large PAR 
defined as: 

W«(i)-fe=£- (2) 

where x^ and a x are respectively the maximum value of 
x n and the root-mean square value of x over all possible 
symbols. 

The real sequence given by Eq. (1) can also be 
obtained by using a complex IFFT of the modulated 
carriers which are organized in the form of the following 
complex sequence Y(n) : 

— — for « = 0, — 
2 2 

Y{n)={X(n) forO<n<y (3) 

0 for— <n£N-\ 
I 2 
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Herein, X{n) is a complex value representing the n*th 
modulated carrier of the multi-carrier symbol to be 
N 

transmitted, — is the number of carriers, and n is an 
2 

integer index. The 1FFT of the complex sequence Y(n) is 
another complex sequence y{n) denoted as: 

y{n)=a{n)+jb{n) (4) 
Herein, a(n)and fc(n)are two real sequences of length N 
and / = . The output y{n)of the IFFT is a complex 
time domain sequence of length N whose real part a{n) 
will be selected for transmission by the real part selector 
REAL. It is further remarked that: 

*( n )=2*(«) = 2ReM*)] (5) 
\fx{n) represents the real IFFT of the multi-carrier 
symbol X(n), i.e. the IFFT of the sub-carrier complex 
conjugate vector built up from X(n). 

When the real part of a sample a(n) of the complex 
time domain multi-carrier symbol y(n) exceeds a certain 
threshold value, the time domain multi-carrier symbol is 
subjected to a transformation. The threshold condition is 
verified by the observation device OBS which thereupon 
controls the transformation unit TRANSF to apply the 
transformation on y{n). Observing the time domain 
multi-carrier symbol y(n) and transforming it is executed 
itcratively. The proposed transformation unit TRANSF is 
able to apply seven different transformations on the 
complex time domain multi-carrier symbol y{n). These 
seven transformations will be discussed further. If none 
of the seven transformations produces a transformed time 
domain multi -carrier symbol whose respective real part 

*i( n )> h( n )> **(*)> 'sM» '*M or h( n ) consist 
of samples with an amplitude below the threshold value, 
a clip event occurs. This may cause harmonics and out- 
of-band radiation on the communication line. Moreover, 
the clipping will result into an increase of the binary error 
rate. 

In a first step, the transformation unit TRANSF 
decomposes the time domain multi -carrier symbol y(n) in 
four partitions each containing information with respect 
to special subsets of carriers of the multi-carrier symbol 
Y(n). Suppose thereto that the frequency domain multi- 
carrier symbol Y{n) would be decomposed into two 
partitions, the first partition Y, carrying the even 
numbered carriers and the second partition Y Q carrying 
the odd numbered carriers. These partitions can 
mathematically be expressed as follows: 



nW=j(l+* 7 "M«) for« = 0,l N-\ (6) 

and 

Y 0 {n)^-e }m \Y{n) for n-O f l,...,N-l (7) 
whereby = Y € {n)+Y 0 {n). 

Y e and Y 0 now can each be decomposed into two further 
partitions: 



\ + e J " 2 ' Y e (n) for n = 0,l N-l (8) 



'A 



f4 



(n) for w = 0,l,. (9) 



whereby Y e (n)= Y Q (n)+Y 2 {n) t and 



\-e*i Y 0 (n) for n = 0, 1 AT — 1 (10) 



Y 0 {n) for n = 0, !,...,#-! (11) 



1\ 



whereby Y u (n) = Y x (n)+ Y) (n) . 

Using the Fourier transform properties, it can be 
verified that the time domain sequences corresponding to 
the previous frequency domain sequences Y e , Y ot Y 0 > 
Y s , Y 2 and Y z are given by: 

y e {n) = a € (n)+jb e {n),wlh (12) 

a r (rt) = I^(n)+a^ff+y) > | for n -0,1 N-\ , and 

and: 

y.W- «.(«)+ A W. with (13) 
fl 0 («)=^ a M~ fl ^ w+ y^ *° r "=0,l,...,W-l ( and 

*«W-^*6i)-*[«+y1i| for n-0 v U. t Ar-l; 



and: 



^0 (») = a 0 i n )+ A ( n ) • W»* ( 1 4 ) 
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and: 



yiW-«lW+^W.wilh (15) 
fl,W-^fl # (n)+^^+ for n=0 > 1,...,AT-l 

and: 

> 2 («) = fl 2 (n)+y'6 2 (n),with (16) 

^W'^eW-^w+jl] for « = 0,U..tf-l 
and: 

> 3 (n)=a 3 (n)+AW,with (17) 

«3W»}[«.W-*.(*+7l] for «=<u....at-i 

friW-^W+^+jV for n=0,l,...,JV-l 

Thus, the transformation unit TRANSF first 
decomposes the time domain multi-carrier symbol y(n) 

in 4 partitions y 0 {n), y\(n), y 2 {n) and y 3 (n). 

Afterwards, the real and imaginary parts of these 
partitions will be used to constitute a transformed time 
domain multi-carrier symbol whose real part is denoted 
by ',M> ' 2 M> h(«)> 'a(»)> h( n )> 'aW or / 7 (/i). 
More specifically, the carriers of a partition must be 
rotated over an angle which is an integer multiple of 



2£ 

N 



radians, the integer multiple being proportional to the 



carrier index n, and will be rotated over a fixed angle 
equal for all carriers of that partition before the partitions 
are re-combined into the transformed time domain multi- 
carrier symbol with real parts /, (n), i = 1, . . . ,7 . 

One of the transformed real time domain symbols 
rj(ft),/s|,...,7 is selected for transmission. It will be 
parallel to serial converted and transferred to the multi- 
carrier receiver. An indication of the applied 
transformation is modulated on the pilot carrier by the 
pilot tone modulator PILOT MOD and added to the 
multi-carrier symbol before transmission towards the 
receiver. Obviously, the applied transformation may be 
reported in another way to the multi-carrier receiver too. 
In the receiver, the transformed time domain multi-carrier 
symbol is again parallelised by the serial to parallel 



converter and converted into a transformed frequency 
domain multi-carrier symbol, respectively 7j(«), T 2 (n), 

r 3 (rt), 7 4 (rt), r s (n), T e (n) or 7 7 («), by the fast fourier 
transfomer FFT, The pilot carrier is monitored by the 
pilot tone demodulator PILOT DMOD, and the 
information demodulated from this pilot carrier is used to 
control the inverse transformation unit INV TRANSF so 
that the exact inverse transformation is applied to the 
transformed frequency domain multi-carrier symbol 
7*/(n),/ = 1,...,7 . Since the inverse transformation is 

executed in the frequency domain, the transformations 
applied by the transformation unit TRANSF in the 
transmitter need to have a simple frequency domain 
equivalent to be usable in practice. When the 
transformation consists of a rotation proportional to the 
carrier index and a fixed rotation equal for all carriers of 
partition, this condition is obviously satisfied. 

In the implementation of the presented technique, 
some specific phasor transformations are used for the 
inverse transformations. Such phasor transformations are 
defined by a phasor transformation vector Z(n). In 
general this phasor transformation vector may be 
expressed as: 

Z(n)=A e Jff ' with n = 0, 1,..., N-\ (18) 
Herein, A is a complex constant with unit magnitude in 
order to not decrease the minimum distance of the 

In 

constellation, and ^ n is any multiple of — proportional 

to the index n, and can be enlarged by a fixed angle g> 0 . 
In other words, <p n can be written as: 

With M being an arbitrary integer and <p 0 being a 
constant angle. Generally, since the time domain multi- 
carrier symbol y(n) is divided into four partitions y 0 (n) y 
y^n), y 2 {n) and ;y 3 (/i), four transformations can be 
defined: Z,(n), Z 2 (n), Z 3 (n)and Z 4 (n). The multi- 
carrier symbol obtained at the output of the carrier 
transformer FFT then can be expressed as: 

T,{n)=Z i {it)-Y 0 ( n )+Z 2 { n )Y l { n )+ 

zUY2(»hzMY 3 (») (20) 

and the time domain multi-carrier symbol transmitted is 
given by: 

f,M- 2 Rc[lFFT(T t («))] (21) 
In the example described below, only two phasor 
transformations Z x (n) and Z 2 (n) are used to perform 
the transformations: 
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z,(«)= 



for n even 



Z^/j^e^-^^ for n odd 
I V2 



(22) 



and: 



r .3^ , . 

/ v J ■ 1—1 

Z 2 \n) = e 4 — =- for /i even 
V2 



(23) 



Z 2 (n)=e J 4 " l^d- for n odd 
I 2 V2 

The seven transformations applied by the 
transformation unit TRANSF on the basis of the phasor 
transformations 2,(/i)and Z 2 (n), and corresponding 
inverse transformations applied by the inverse 
transformation unit INV TRANSF are listed hereafter. 

I. According to the first transformation, the real and 
imaginary parts of the four partitions j>oM> >*iM» 
y 2 {n) and>> 3 (/i)are combined by the transformation unit 
TRANSF into the following transformed time domain 
multi-carrier symbol: 



, («)= 2{a 0 (/.)+ a 3 (/0)+V^ 2 (* + jl+ 



If the corresponding transformed frequency domain 
multi-carrier symbol, i.e. the multi-carrier symbol 
obtained at the output of the fourier transformer FFT in 
the receiver RX when t x (n) is supplied thereto, is 
denoted by T x (n) , the original frequency domain multi- 
carrier symbol can be recovered from T { (n) by the 
inverse transformation unit INV TRANSF if the 
following inverse transformation is applied: 

^2r,(/i) for «=o,y 

T x (n ) for n € {0 mod 4, 3 mod 4} and 

2 

*<»HHW for „ e { 2mod 4}and2 S *^-l< 25 > 
Z x \n) 2 

^^forne{lmod4}and !<*<— 1 
Z 2 \n) 2 



r 



for— <n<N 
2 



Herein, the notation {k mod/} stands for the set 
{*, k + /, k + 21, k + 3/,...} , e.g. {0 mod 4} stands for the set 
{0,4,8,12,...}. 

II. According to the second transformation, the four 
partitions y$(n\> v,(«), y 2 (n) and ^(/i) are combined 
by the transformation unit TRANSF into the following 
transformed time domain multi-carrier symbol: 



r 2 (/ J )=2(a 0 (^Ka l ( W ))+^ 2 ^ + ^ > |+ 



(26) 



If the corresponding transformed frequency domain 
multi-carrier symbol is denoted by T 2 (n), the original 
frequency domain multi-carrier symbol can be recovered 
from T 2 (n) by the inverse transformation unit if the 
following inverse transformation is applied: 

^2T 2 {n) for n = 0,y 

7* 2 (n ) for n e {0 mod 4, 1 mod 4}and 

0<Sn£ — -1 
2 

(24 ) * W - ILto for n 6 { 2 mod 4 }and 2 S n* 1 < 27 > 



^4 for *e{3rnod4}and3</j<— -1 
Z 2 (n) 2 



,0 



N 

for — < n< N 
2 



III. 



L 

In the third transformation, the four 
partitions y 0 (n) , y x (n) , y 2 {n) and y 3 {n) are combined by 
the transformation unit TRANSF into the following 
transformed time domain multi -carrier symbol: 



3 (n)=2^(n)+>/2^d l ^+jV6 l ^ + ^lH 



(28) 



If the corresponding transformed frequency symbol is 
denoted by T 3 (n) 9 the original frequency domain multi- 
carrier symbol can be recovered fromr 3 (/i) by the 
inverse transformation unit if the following inverse 
transformation is applied: 
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2T 3 {n) for /i = 0,y 

f 3 (n) for n even 

jf^^IlW for w€ {lmod4}andl£«£ — 1 (29) 
Z,(/i) 2 

for /!<= {3mod4}and3£n <S — -1 
2 2 («) 2 

0 for— 

1 2 

IV. According to the fourth transformation, the four 
partitions yo{n) t y l (n) i y 2 {n)^nd y 3 {n) are combined 
into the following transformed time domain multi-carrier 
symbol: 



(30) 



At the receiver, the original frequency domain multi- 
carrier symbol can be recovered from the corresponding 
transformed frequency multi-carrier symbol if the 
following inverse transformation is applied: 

^2T A {n) for « = 0,y 

T 4 {n) for #te {0mod4,3mod4}and 

*MHl40 for ^ 

Z 2 [n) 2 
%H for n e {l mod 4}and 1 S n Z 1 

.0 for — <n<N 

I 2 
V. In the fifth transformation, the four partitions 
yo i n )*y\{ n )*yi{ n )* n & y*{ n ) are combined into the 
following transformed time domain multi-carrier symbol: 



'5 W = 2(a 0 („)+ (*))+ V5ja 2 (* +^) + 



(32) 



The original DMT symbol A"(w)can be recovered 
from the corresponding frequency symbol ^(njby 
applying the following transformation: 



^2T s (n) fom = 0,y 
T 5 {n) for n € {0 mod 4, 1 mod 4}and 



0£n£— -1 
2 



*(«H2iW for „ e {2mod4}and2 < w < ^_,(33) 

for ne {3mod4}and3 <n £ — -1 
Zj^/i) 2 



for — </i<# 
2 



VI. In the sixth transformation, the four partitions 
yoi n )>y) M» J^M*™* .VjM arc combined into the 
following transformed time domain multi-carrier symbol: 



6 (rt) = 2a,(«)+^,^+^V&,^+^lH 



(34) 



The original DMT symbol A^njcan be recovered 
from the corresponding frequency symbol T 6 (n)by 
applying the following transformation: 

r 

2T 6 {n) for/i = 0,y 
T 6 (n) for n even 



xM^i^rl for #ie{linod4}andl (35) 
Z 2 \n) 2 

for /ie{3mod4}and3:S/i£— -1 



for —<n<N 



VII. The seventh transformation is obtained by a 
simple modification of the previous transformation, the 
corresponding time domain multi-carrier symbol is given 
by: 



(36) 



The original frequency domain DMT symbol X(n)c*n 
be recovered from the output of the fourier transformer 
FFT in the receiver, denoted 7* 7 (/i), corresponding to the 
input r 7 (n)if the following transformation is applied: 
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f 



2r 7 (/i) for/i = 0,y 
7*7 (n) for w even 

'TP 

-^7 for ne {3mod4}and3 < n <>— -1 

for — <n<W 
2 



X(n) = C Tl ft for /ie{lmod4}andl£/i<— -1 (37) 
Z 2 (/i) 2 



r 



3. Complexity and gain 

At the transmitter, one needs to compute a preliminary 
output to detect if the algorithm must be applied or not, 
this extra information is encoded into a pilot and send to 
the receiver. If the proposed algorithm was applied, the 
complexity at the transmitter measured in terms of extra 
number of additions and multiplications is limited by 
(5+3*3+2*3)N real additions for T u T 2i T 4 and T s and 
(5+3*4+2*2)N real additions for r 3 ,r 6 and T 7 and 2N 

real multiplications for each transformation. At the 
receiver, the original DMT symbol can be recovered by 

performing maximum ^ complex multiplications at the 

output of FFT transformer. 

One can mention that the seven transformations 
described above are an example corresponding to the 
most complex case. There is still room to improve this 
technique by decreasing the implementation complexity 
as will be shown in the next section. 

Another remark is that the proposed algorithm is not 
restricted to ADSL systems. It can be used for all systems 
wherein a multi -carrier modulation is applied with QAM- 
constellations of different size used in each carrier. 
Figure 1 shows the clip probability versus the allowed 
PAR for an ADSL implementation using N-512 and 16- 
QAM in each carrier. The curves denoted by R4 and R8 
correspond to Orckit method [7] when applying the 
original IFFTs and three or seven transformations 
respectively. Similarly, the curves denoted by S4 and S8 
correspond to applying the original IFFTs and 
transformations 1, 11, 111 or the complete transformations 
respectively. The maximum gain in PAR introduced by 
the proposed technique is about 3.7 dB instead of 2.7 dB 
given in [7]. 

The simulation results show that a given maximum 
clipping probability can be satisfied whose power level is 
at least 1.5 dB below the clipping power level obtained 
by the technique described in [7] corresponding to the 
same maximum clipping probability. 




10 11 12 13 14 15 



Figure 1. Clipping probability 
4. Reducing the implementation complexity 

In this section we propose some modifications to 
reduce the complexity of the previous algorithm. As can 
be observed, in each transformation described above, 

there are two scale factors 2 and in order to 

eliminate the factor V2 , two possible solutions can be 
proposed: 

The first one is to choose a suitable phasor 
transformation Z,(/i)and Z 2 {n) as for example: 



I j-n 



for n even 



* 2 for n odd 



and: 



r a. 



for n even 



(38) 



(39) 



for n odd 

In this case, the number of multiplications at the 
transmitter will be reduced by 50% at the expense of an 
increase of the number of additions by an amount 2N per 
transformation. 

The second one is to rotate all the four partitions 
Yo( n )> Y\( n )> l^Mand J^Mby a specific phasor 
transformation Z,(/t) as defined in Eq. (18) and then a 
new frequency sequence can be given by: 
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r < («)= z a (4Y 0 (n) + zMYd"hzMY2(») Mm 

The complexity of this solution is exactly the same as 
in the previous case with a difference that the entire 

symbol now is scaled by a factor instead of 2. 
The simulation results show that the improvements in 
PAR reduction given by this solution are almost the same 
given in figure 1 . 

5. Conclusion 

A new method for clip probability reduction is 
proposed by reducing the PAR of multi-carrier signal. 
This method is based on the shift property of an FFT in 
the time domain. The new time domain signals can be 
easily computed at the transmitter and the receiver can 
efficiently recover the original frequency symbol from 
the received frequency symbol. A gain of more than 3.7 
dB can be achieved with low complexity and redundancy. 
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